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ABSTRACT 

The Hubble Ultra Deep Field (HUDF) contains a significant number of i?-, V- and i'-band dropout 
objects, many of which were recently confirmed to be young star-forming galaxies at 0~4— 6. These 
galaxies are too faint individually to accurately measure their radial surface brightness profiles. Their 
average light profiles are potentially of great interest, since they may contain clues to the time since 
the onset of significant galaxy assembly. We separately co-add V, i'- and z'-band HUDF images of sets 
of 2;~4, 5 and 6 objects, pre-selected to have nearly identical compact sizes and the roundest shapes. 
From these stacked images, we are able to study the average(d) radial structure of these objects at 
much higher signal-to-noise ratio than possible for an individual faint object. Here we explore the 
reliability and usefulness of a stacking technique of compact objects at z~4 — 6 in the HUDF. Our 
results are: (1) image stacking provides reliable and reproducible average surface brightness profiles; 
(2) the shape of the average surface brightness profile shows that even the faintest z~4 — 6 objects 
are resolved] and (3) if late-type galaxies dominate the population of galaxies at z~4— 6, as previous 
HST studies have shown for z<4, then limits to dynamical age estimates for these galaxies from their 
profile shapes are comparable with the SED ages obtained from the broadband colors. We also present 
accurate measurements of the sky-background in the HUDF and its associated Icr uncertainties. 
Subject headings: galaxies: high-redshift — galaxies: structure — galaxies: formation 



1. INTRODUCTION 

In the last decade, ground and space based observa- 
tions of high redshift galaxies have begun to outline the 
process of galaxy assembly. The details of that process 
at high redshifts, however, remain poorly constrained. 
There is increasing support for the model of galaxy for- 
mation, in which the most massive galaxies assemble ear- 
lier than their less massive counterparts (e.g. Cowie et al. 
1996; Guzman et al. 1997; Kodama et al. 2004; McCarthy 
etal. 2004). A detailed analysis of the 'fossil record' 
of the current stellar populations in nearby galaxies se- 
lected from the Sloan Digital Sky Survey (SDSS; York 
et al. 2000) provides strong evidence for this downsizing 
picture (Heavens etal. 2004; Panter etal. 2007). The in- 
creasing number of luminous galaxies spectroscopically 
confirmed to be at z^Q.h (e.g. Hu etal. 2002; Kodaira 
etal. 2003; Kurk etal. 2004; Rhoads etal. 2004; Stern 
etal. 2005; Taniguchi etal. 2005), or <0.9 Gyr after the 
Big Bang, also supports this general picture. In an al- 
ternate hierarchical scenario, arguments have been made 
that significant number of low luminosity dwarf galaxies 
were present at these times, and were the main contribu- 
tor to finish the process of reionization of the intergalac- 
tic medium (Yan & Windhorst 2004a, b). However, there 
is presently little information on the dynamical struc- 
ture of these or other galaxies at z ~ 6. It is not clear 
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whether these objects represent isolated disk systems, or 
collapsing spheroids, mergers or other dynamically young 
objects. 

Ravindranath et al. (2006) used deep, multi- 
wavelength images obtained with the Hubble Space 
Telescope [HST] Advanced Camera for Surveys (ACS) 
as part of the Great Observatories Origins Deep Survey 
(GOODS) to analyze 2-D surface brightness distribu- 
tions of the brightest Lyman Break Galaxies (LBGs) 
at 2.5 < z < 5. They distinguish various morphologies 
based on the Sersic index n, which measures the shape 
of the azimuthally averaged surface brightness profile 
(where n=\ for exponential disks and n=4 for a de 
Vaucouleurs law). Ravindranath etal. (2006) find that 
40% of the LBGs have light profiles close to exponential, 
as seen for disk galaxies, and only ~30% have high n, 
as seen in nearby spheroids. They also find a significant 
fraction (^30%) of galaxies with light profiles shallower 
than exponential, which appear to have multiple cores 
or disturbed morphologies, suggestive of close pairs or 
on-going galaxy mergers. Distinction between these 
possible morphologies and, therefore, a better estimate 
of the formation redshifts of the systems observed at 
2~4— 6 in particular, is important for testing the galaxy 
assembly picture, and for the refinement of galaxy 
formation models. 

One possible technique involves the radial surface 
brightness profiles of the most massive objects — those 
that will likely evolve to become the massive elliptical 
galaxies, which we see in place at redshifts z < 2 (Driver 
etal. 1998; van Dokkum etal. 2003, 2004). This can be 
analytically understood in the context of the Lynden-Bell 
(1967) relaxation formalism and the numerical galaxy 
formation simulations of van Albada (1982), which de- 
scribe collisionless collapse and violent relaxation as the 
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formation mechanism for elliptical galaxies. As the time- 
scale for relaxation is shorter in the inner than in the 
outer parts of a galaxy, convergence toward a r^/"-profile 
will proceed from the inside to progressively larger radii 
at later times. Moreover, Kormendy (1977) has shown 
that tidal perturbations due to neighbors can cause the 
radial surface brightness profile to deviate from a pure 
de Vaucouleurs profile in the outer parts of a galaxy. This 
implies that the radius where surface brightness profiles 
start to deviate significantly from an r^l"^ profile might 
serve as a "virial clock''' that traces the time since the on- 
set of the last major merger, accretion events or global 
starburst in these objects. 

Image stacking methods have been used extensively 
on X-ray (Nandra etal. 2002; Brandt etal. 2001) and 
radio (Georgakakis et al. 2003; White et al. 2007) data to 
study the mean properties (e.g. flux, luminosity) of well- 
defined samples of sources that are otherwise too faint to 
be detected individually. Pascarelle et al. (1996) applied 
such a stacking method to a large number of optically 
very faint, compact objects at 2; = 2.39 to trace their 
"average" structure. This approach was also applied by 
Zibetti et al. (2004) to detect the presence of faint stellar 
halos around disk galaxies selected from the SDSS. An 
attempt to apply this technique to high redshift galaxies 
in the Hubble Deep Field (HDF; Williams et al. 1996) was 
not conclusive (H. Ferguson; private communication) due 
to the poorer spatial sampling and shallower depth of the 
HDF compared to the Hubble Ultra Deep Field (HUDF; 
Beckwith etal. 2006). 

In this paper, we use the exceptional depth and fine 
spatial sampling of the HUDF to study the potential of 
this image stacking technique, and will estimate limits 
to dynamical ages of faint, young galaxies at z ~ 4 — 6. 
The HUDF reaches ^^1.5 mag deeper than the equiva- 
lent HDF exposure in the i'-band, and has better spatial 
sampling than the HDF. The HUDF depth also allows 
us to characterize the sky background very accurately, 
which is critical for successfully using a stacking method 
to measure the mean surface-brightness profiles for these 
faint young galaxies. 

This paper is organized as follows: In § 2 we summarize 
the HUDF observations, and in § 3 we discuss the selec- 
tion of our z~4, 5 and 6 samples. In § 4 we describe our 
data analysis, which includes accurately measuring the 
Ict sky-subtraction error, the image stacking method to 
generate mean surface-brightness profiles, and our test of 
its reliability. In § 5 we present and discuss our results in 
terms of the average surface-brightness profiles of 2 ~ 4^6 
galaxies, and in § 6 we conclude with a summary of our 
results. 

Throughout this paper we refer to the HST/ AGS 
F435W, F606W, F775W, and F850LP filters as the B-, 
V-, i'-, and z'-bands, respectively. We assume a Wilkin- 
son Microwave Anisotropy Probe (WMAP) cosmology of 
n„=0.24, f2A=0.76andFo=73kms~^ Mpc~^, in accord 
with the most recent 3-year WMAP results of Spergel 
etal. (2007). This implies a current age for the Universe 
of 13.65 Gyr. All magnitudes are given in the AB system 
(Oke & Gunn 1983). 

2. OBSERVATIONS 

The HUDF contains >100 objects that are i'-band 
dropouts, making them candidates for galaxies at 2; ~ 6 




z' (AB mag) 

Fig. 1. — The HUDF number counts for all 2: ~ 4, 5, 6 objects 
before the sub-selection of compact isolated 2 ~ 4, 5,6 objects was 
made. The vertical dotted line shows the magnitude to which the 
number counts of all these redshifts are complete. The area of the 
HUDF is 3.15x10-3 deg^. 

(Bouwens et al. 2004, 2006; Bunker & Stanway 2004; Yan 
& Windhorst 2004b). Similarly, there are larger num- 
bers of objects in the HUDF that are B-band dropouts 

(415 in total) or U-band dropouts (265 in total), and 
are candidates for galaxies at z ^ 4 and 2^5, re- 
spectively. Beckwith et al. (2006) and Bouwens et al. 
(2007) find similar number of B- and U-band dropouts 
in the HUDF. A significant fraction of these objects to 
AB<27 mag have recently been spectroscopically con- 
firmed to have redshifts z ~ 4 — 6 through the detection 
of Lya emission or identifying their Lyman break (Mal- 
hotra etal. 2005; Dow-Hygelund etal. 2007). We discuss 
our detailed drop-out selection criteria below. Despite 
the depth (AB<29.5 mag) of the HUDF images, how- 
ever, these objects appear very faint and with little, if 
any, discernible structural detail. Visual inspection of 
all these objects shows their morphologies to divide into 
four broad categories: symmetric, compact, elongated, 
and amorphous. 

3. SAMPLE SELECTION 

We constrtict three separate catalogs for these z ~ 
4, 5, 6 galaxy candidates, selecting only the isolated, 
compact and symmetric galaxies. We exclude objects 
with obvious nearby neighbors, to avoid a bias due to 
dynamically disturbed objects and complications due 
to chance superpositions. Figure 1 demonstrates that 
our completeness limit for 2; ~ 4 and 2; ~ 5 ob- 
jects is AB<29.3 mag, and for z ~ 6 objects it is 
AB<29.0 mag. Therefore, all three catalogs are com- 
plete to AB<29.0 mag, which is equivalent to at least a 
10(7 detection for objects that are nearly point sources. 
For each object in our 2; ~ 4, 5, 6 samples, we extracted 
51x51 pixel postage stamps (which at 0'.'03 pix"'^ span 
1'.'53 on a side) from the HUDF V , i' and z'-band images, 
respectively. Each postage stamp was extracted from the 
full HUDF, such that the centroid of an object (usually 
coincident with the brightest pixel) was at the center of 
that stamp. 
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3.1. The 2~4 andzc^h Objects (B-, V-band 
dropouts) 

We used the z'-band selected BVi'z' HUDF catalog 
(Beckwith et al. 2006) to select the z ~ 4 and z ~ 5 ob- 
jects. With the HyperZ code (Bolzonella et al. 2000), 
we computed photometric redshift estimates, using the 
magnitudes and associated uncertainties tabulated in the 
HUDF catalog. All objects with 3.5< -2phot ^4.5 were as- 
signed to the bin of 2 ~ 4 candidates, and all objects with 
4.5< Zphot <5.5 to the bin of z~5 candidates. 

We then applied color criteria, similar to those adopted 
by Giavalisco etal. (2004), to select the ~ 4) and 
V{z ~ 5) dropout samples. For B-band dropouts, we 
require: 

( {B -V) >1.2 + lAx (y - z') mag 
<^ and iB-V)> 1.2 mag 
i and {V - z') < 1.2 mag 

For F-band dropouts, the following color selection was 
applied: 



(V-i') > 1.5 + 0.9 X (j'- 
and {V - i') > 1.2 mag 
and {i' — z') < 1.3 mag 



z') or {V - i') > 2.0 mag 



We note, that only objects satisfying both color and pho- 
tometric redshift criteria were selected in our samples. 
Vanzella etal. (2006) using VLT/F0RS2 observed ~100 
B-, V- and i'-band dropout objects in the Chandra Deep 
Field South (CDFS) selected based on above mentioned 
color criteria (Ciavalisco et al. 2004). They have spcctro- 
scopically confirmed >90% of their high redshift galaxy 
candidates. Therefore, we expect only a small number 
(<10%) of contaminants in our sample of dropouts. One 
or two objects in our final sample could be such contam- 
inants, but because we have 3 different realizations of 
10 objects (3x10), each showing similar profiles, they do 
not appear to affect our results. 

The z~4 sample has 415 objects, while the z~5 sam- 
ple has 265 objects. In Figure 2ab, we show the distribu- 
tion of the FWHM and ellipticity, e = (1 — 6/a), measured 
in each of the two samples using SExtractor (Bertin 
& Arnouts 1996). We further constrained our samples 
by imposing limits on compactness and on roundness of 
FWHM < 0'.'3 and e < 0.3. Again, this is to minimize the 
probability that the zc^A—b candidates are significantly 
dynamically disturbed, and to maximize the probability 
of selecting physically similar objects. Our goal is to find 
the visibly most symmetric, least disturbed systems for 
the current study. This sub-selection leaves 204 objects 
in the 2^4 sample and 102 objects in the 2:~5 sample. 
Most of these objects are faint, and are only a few pix- 
els across in size, and, hence, have larger uncertainties 
in their measurements of FWHM and ellipticity. There- 
fore, we also checked our objects visually to eliminate any 
possibility of our selected objects being contaminated by 
unrelated nearby objects, being clearly extended, or ob- 
jects with complex morphologies. 

3.2. 2;~6 Objects (i' -band dropouts) 

Yan & Windhorst (2004b) found 108 possible 5.5< 
z<6.5 candidates in the HUDF to mABiz85o)='^Q-0 mag. 
Bunker & Stanway (2004) independently found the 
brightest 54 of these 108 z ~ 6 candidates to 
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Fig. 2. — Ellipticity, (1 - b/a), versus object FWHM, for all 
2~4 (a), 2~5 (b) and z~6 (c) objects selected in the HUDF. 
Measurements were performed in j'-band for ^ ~ 4 and z ~ 5 ob- 
jects, while we used the 2:'-band for z ~ 6 objects. The FWHM 
of a stellar image/PSF is ~3 pixels or (y.'09, indicated by the left- 
most hatched area in each panel. Objects within the shaded area 
meet our additional selection criteria on roundness (e < 0.3) and 
compactness (FWHM < (/.'S or 10 pixels). 



AB=28.5 mag. Similarly, deep HST/ ACS grism spectra 
of the HUDF i'-band dropouts confirm >90% of these 
objects at AB<27.5 mag to be at z~6 (Malhotra etal. 
2005; Hathi etal. 2007). Using the catalog of Yan & 
Windhorst (2004b), we extracted 108 postage stamps, 
each 51x51 pixels in size, from the HUDF z'-band im- 
age. 

Like for the z ~ 4 and z ~ 5 objects, for each z ~ 6 
object we measured its z'-band FWHM and ellipticity 
using SExtractor. Figure 2c shows the measured ellip- 
ticity versus FWHM for all 108 z ~ 6 candidates. A 
smaller sample of 67 objects satisfies our constraints on 
the FWHM and ellipticity. Further visual inspection, to 
make sure that our sample has only isolated, compact 
and rotmd objects, leaves 30 objects in our z~6 sample. 
We therefore imposed a sample size of 30 objects also on 
the two lower redshift bins after visual inspection. 

The results in this paper are therefore based 
on approximately (30/415)~7%, (30/265)^11%, and 
(30/108)~28% of the total z ~ 4,5 and 6 galaxy pop- 
ulations. 

4. RESULTS 
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4.1. The HUDF Sky Surface- Brightness Level and its 
rms Variation 

For the present work, it is critical that we accurately 
characterize the sky-background, and correctly propa- 
gate the true la errors due to the subtraction of this 
sky-background. In the following, we will pursue two 
complimentary approaches to determine the sky surface- 
brightness, and compare the results. Here, we discuss 
the 2;'-band measurements in detail. 

We first measured the sky-background in each of 
the 415 z ~ 4 object stamps ('local' sky measure- 
ments). The Interactive Data Language (IDL^) pro- 
cedure SKY/MMM.pro^ was used to measure the sky- 
background. This procedure is adapted from the 
DAOPHDT (Stetson 1987) routine of the same name and 
works as follows. First, the average and sigma arc ob- 
tained from the sky pixels. Second, these values are used 
to eliminate outliers with a low probability. Third, the 
values are then recomputed and the process is repeated 
up to 20 iterations. If there is a contamination due to an 
object, then the contamination is estimated by compar- 
ing the mean and median of the remaining sky pixels to 
get the true sky value. The output of this procedure is 
the modal sky-level in the image. 

Figure 3c shows a histogram of the z'-band modal sky 
values obtained from all 415 object stamps extracted 
from the drizzled HUDF images. The la uncertainty 
in the sky, crgky, determined from a Gaussian fit to 
the histogram, is 2.19 x 10~^ electrons sec~^ in the z'- 
band. The sky-background level within the HUDF was 
obtained from the original flat-fielded ACS images, be- 
cause the final co-added HUDF data products are sky- 
subtracted. The header parameters MDRIZSKY and 
EXPTIME were used to obtain the actually observed 
sky- value. MDRIZSKY is the sky value in electrons (e~) 
computed by the MultiDrizzle code (Kockcmoer et al. 
2002), while EXPTIME is the total exposure time for 
the image in seconds, so that the average sky- value in the 
HUDF has the units of e~ sec~^. Figure 4d shows the 
histogram of the sky- values obtained from 288 HUDF z'- 
band flat-fielded exposures. The average value of the sky 
background, Isky, is 0.02051 e~ sec"^ pix^^. That sky- 
value is measured from the flat-fielded individual ACS 
images with pixel sizes of 0'.'05 pix~^ and hence, in the 
following calculations, the average sky- value is multiplied 
by a factor of (0.030/0.05)2=0.602 to obtain the corre- 
sponding average sky- value for the HUDF drizzled pixel 
size of O'.'OSO pix~^. Using these values, we estimate the 
relative rms random sky-subtraction error as follows: 



2.19 X 10-5 



I, 



sky 



2.05 X 10-2 • 0.602 



= 2.97 X 10" 



The measured average sky background level can then 
be expressed as the 2:'-band sky surface brightness as 
follows: 



H^, = 24.862 - 2.5 • log 



0.0205 • 0.60^ 
0.0302 



= 22.577 ± 0.003 mag arcsec-^ 

^ IDL Website http://www.ittvis.com/index.asp 
^ Part of the IDL Astronomy User's Library, 
http://idlastro.gsfc.nasa.gov/homepage.html 
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Fig. 3. — Distribution of the modal sky background level used to 
estimate the Ict uncertainty in that level, as measured in the 415 
2~4 object stamps extracted from the drizzled HUDF images (a) 
for V-band, (b) for i'-band and (c) for ^'-band. The mean 
and the sigma (cr) of the best-fit Gaussian to these distributions 
are also shown in each panel. 
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M=0.02051 
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Fig. 4. — The actual sky values measured using header param- 
eters MDRIZSKY and EXPTIME from fiat-fielded HUDF expo- 
sures, (a) for B-baiid using 112 exposures, (b) for V-band using 
112 exposures, (c) for i'-band using 288 exposures, (d) for z'- 
band using 288 exposures. The mean (fi) and the sigma (cr) of the 
best-fit Gaussian to these distributions are shown in each panel. 
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where 24.862 is the ACS/WFC z'-band AB zero-point, 
and O'.'OSO pixcl^^ is the drizzled pixel scale. This is 
consistent with the values obtained by extrapolating the 
on-orbit BVI sky surface brightness of Windhorst et al. 
(1994, 1998) to z', with the sky-background estimates 
from the ACS Instrument Handbook (Gonzaga et al. 
2005), and with the colors obtained by convolving the 
filter transmission curves with the solar spectrum. Ta- 
ble 1 gives the measured electron detection rate, sur- 
face brightness and colors of the sky background with 
their corresponding errors for the HUDF BVi'z' bands 
as calculated from Figure 3 and Figure 4. The contri- 
bution of the zodiacal background dominates the total 
sky-background, which we find to be only ^10% redder 
in {V-i') and {i'-z') than the Sun. The ^'-band sur- 
face brightness corresponding to the la sky-subtraction 
uncertainty is therefore: 

- 2.5 • log(S,s,ran) = 22.577 - 2.5 • log(2.97 x IQ-^) 

= 28.895 mag arcsec"^ 

Next, we measure the sky-background from 415 'blank' 
sky stamps (51x51 pixel) distributed throughout the 
HUDF ('global' sky measurements). We measure the 
sky background using the same IDL algorithm as used 
above. 

Figure 5c shows the histogram of the measured z'-band 
modal sky values. A Gaussian distribution was fit to this 
histogram, giving a sky-sigma of 2.00 x 10~^ e~ sec~^. 
The average value of the sky remains 0.02051 sec~^ 
(Figure 4d). Using these values, we can estimate a rela- 
tive rms systematic sky-subtraction error as follows: 



2.00 X 10-5 



I 



sky 



2.05 X 10-2 • 0.602 



= 2.71 X 10-^ 



Since the z'-band sky surface brightness remains 22.577 
mag arcsec-2, this gives us for the surface brightness 
corresponding to the la sky subtraction uncertainty: 

/i^/ - 2.5 • log(Sss,sys) = 22.577 - 2.5 • log(2.71 x 10-^) 
= 28.995 mag arcsec-^ 

From these two complementary approaches, we can 
conclude that all surface brightness measurements be- 
come unreliable for surface-brightness levels fainter than 
28.95±0.05 mag arcsec-^ in the z'-band. We have also 
experimented with slightly larger cutouts (75x75 pixels 
instead of 51x51 pixels) to estimate the sky-subtraction 
error. We find that with the larger cutouts, the surface 
brightness corresponding to the la sky-subtraction er- 
ror is ~0. 1-0.2 mag arcsec-^ fainter. For larger cutouts 
we expect this surface brightness to be ~0.4 mag fainter 
but we find about 0.1-0.2 mag fainter. This might be be- 
cause of residual systematic errors in the HUDF images. 
Therefore, we are at the limit of accurately measuring 
this surface brightness and hence, we will here quote the 
conservative brighter limit of the surface brightness cor- 
responding to this la sky-subtraction error. Expected 
contributions to this surface brightness due to uncertain- 
ties in the bias level determinations, which correspond to 
~0.001 counts sec-^ for typical HUDF exposures (A. M. 
Koekemoer; private communication), are less than 1%. 

Figure 5 clearly shows that the distribution of the 
modal sky-values is not as symmetric around zero as 
in Figure 3, and hence, the use of a 'global' sky value 
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Fig. 5. — Distribution of the modal sky background level used 
to estimate the la uncertainty in that level, as measured in 415 
'blank' 51x51 pixel sky stamps extracted from the drizzled HUDF 
images (a) for V-hand, (b) for t'-band and (c) for z'-band. The 
mean (/i) and the sigma (cr) of to the best-fit Gaussian to these 
distributions are shown in each panel. 

for the HUDF is not as reliable as 'local' sky measure- 
ments. Therefore, for the surface brightness profiles and 
the following discussion, we will adopt the local la ran- 
dom sky-subtraction error for all objects in our study. 

The average modal sky values and their la errors in the 
V- and j'-bands were calculated in exactly the same way 
as for the z'-band, as shown in Figure 3, 4 and 5. The re- 
sulting BVi'z' sky values and the sky surface-brightness 
levels are all given in Table 1. 

4.2. Composite Images and Surface Brightness Profiles 

For each redshift bin (z~4, 5,6), we generated three 
"stacked" composite images from subsets of 10 postage 
stamps that were selected as follows. After placing all 30 
image stamps per redshift bin into a30x (51x51) pixel 
IDL array. 10 stamps were randomly drawn without se- 
lecting any object more than once. An output image 
was generated, in which the values at each pixel are the 
average of the corresponding pixels in the 10 selected in- 
put stamps. From the remaining 20 stamps, we again 
randomly select 10, from which we generated a second 
composite image, after which the final 10 images were 
averaged into the third composite image. The three com- 
posite images per redshift bin are therefore independent 
of each other. In none of our realizations did we pro- 
duce composite images that were essentially unresolved. 
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TABLE 1 

Measured sky values in BVi'z' (filters) for the HUDF 



HUDF 


Number of 


Mean Sky Value" 


Sky SB'^ 




Sky Color'= 


1(7 Sky-Subtraction 


Filter 


Exposures 


(e~/s) and rms error*" 


(AB mag arcsec" ^) 




(AB mag) 


error (AB mag arcsec"^) 


B 


112 


0.015909 ± 0.000065 


23.664 ± 0.003 


(B 
(V 


- V)sky =0.800 


29.85 ± 0.05 


V 


112 


0.070276 ± 0.000297 


22.864 ± 0.002 


- i')sky=0.222 


30.15 ± 0.15 


i' 


288 


0.040075 ± 0.000088 


22.642 ± 0.002 


ii' 


- 2').ky=0.065 


29.77 ± 0.20 


z' 


288 


0.020511 ± 0.000047 


22.577 ± 0.003 


{V 


- 2')sky=0.287 


28.95 ± 0.05 



^ From Figure 4 

^ Error is standard deviation of the mean {cr/^N) 

^ Sky surface briglitness values and colors are consistent with the solar colors in AB mag of (y-i')=0.19, (V-z')— 0.21 and 
{i'-z')=0.01 [except for bluest color (B-V)], and is dominated by the zodiacal background. 



Even the faintest 2:~4 — 6 galaxies are clearly resolved. 
The 2; ~ 4, 5, 6 objects used to generate the composite 
images have an apparent magnitude range of approxi- 
mately 27.5±1.0 AB mag. Because the magnitude range 
is relatively small and the S/N per pixel is low even in 
their central pixel, we have given all objects equal weight. 
To test whether this range in magnitude will affect our 
stacks and hence, our profiles, we created 3 stacks de- 
pending on the apparent magnitude, i.e. one stack of 
the 10 brightest objects in the sample, a second stack of 
the 10 next brightest objects in the sample and a third 
stack of the 10 faintest objects in the sample. This is 
summarized in Figure 6d. We found that the profiles 
were very similar except that the profiles of the fainter 
stacks fall-off more quickly at larger radius compared to 
the profile of the brightest stack, but the inner profile 
and the deviation in the profiles are clearly visible in 
all 3 stacks. Therefore, we conclude that for our range 
in apparent magnitudes, our stacks/profiles are not af- 
fected. Perhaps most surprisingly. Figure 6d shows that 
re value of all 3 flux ranges (-26.0-27.0, -27.0-28.0 & 
—28.0-29.0 mag) are all about the same over —3-4 mag 
in flux, so the primary parameter that distinguishes the 
brighter from the fainter 2 c± 6 dropouts is their central 
surface brightness (which thus also varies by —3-4 mag). 

We used the IRAF'^ procedure ELLIPSE to fit sur- 
face brightness profiles shown in Figure 6 to each of the 
three independent composite images per redshift bin. We 
also computed a mean surface-brightness profile from 
the three composite surface brightness profiles gener- 
ated from the three independent composite images for 
each redshift bin. Figure 7 shows composite images for 
2 — 4, 5, 6 objects. Here each composite image is a stack 
of 30 objects. Figure 8 shows the average surface bright- 
ness profiles for each of the redshift intervals 2; ~ 4, 5, 6. 
The thin solid curves in Figure 6 and the dot-dash curves 
in Figure 8 represent the observed ACS V , i' and z'- 
band Point Spread Functions (PSFs), while the horizon- 
tal dashed lines indicate the surface brightness level cor- 
responding to the 1(7 sky-subtraction error in each of the 
HUDF images as discussed in § 4.1. It is important to 
note that we scaled the ACS PSFs to match the surface 
brightness of the central data point in our mean surface- 
brightness profile, to determine how extended the mean 
surface-brightness profile is with respect to the PSFs. 

''' IRAF (http://iraf.net) is distributed by the National Optical 
Astronomy Observatories, which are operated by the Association 
of Universities for Research in Astronomy, Inc. , under cooperative 
agreement with the National Science Foundation. 
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Fig. 6. — Composite surface brightness profiles for three indepen- 
dent sets of 10 objects at (a) z~4, (b) z~5 and (c) ^~6, respec- 
tively. The thin solid curve represents the ACS V, i' and z'-band 
PSFs, respectively, while the horizontal dashed line indicates the 
surface brightness level corresponding to the 1<t sky-subtraction 
error in the HUDF images. The vertical dotted line marks the 
radius at which the profile starts to deviate significantly from the 
extrapolation of the inner r^/" profile observed at smaller radii. 
Note that at z ~ 6, this deviation is still well above the red z'- 
band PSF halo at r>0'.'30. The panel (d) shows 3 z~6 composite 
profiles (each with a set of 10 objects) divided by apparent mag- 
nitudes. The brightest composite profile (dotted) has an average 
z'-band magnitude of ~26.8 mag. The next brightest composite 
profile (short dash) has an average z'-band magnitude of ~27.9 
mag, and the faintest composite profile (dot-dash) has an average 
z'-band magnitude of ~28.9 mag. 

In Figure 8, we fitted all possible combinations of the 
Sersic profiles (convolved with the ACS PSF) to the ob- 
served profiles and using minimization, found the best 
fits for galaxies at 2; ~ 4, 5, 6. The best fit Sersic index 
(n) for all three profiles (z — 4, 5, 6) is n < 2, mean- 
ing these galaxies follow mostly exponential disk-type 
profiles in their central regions. We find that the ob- 
served profiles start to deviate from the best-fit profiles 
at r >0'.'27, somewhat depending on the redshift. Prom 
Figure 8, we also sec that in each of (z — 4), i' (z — 5) 
and z' (z — 6), the PSF declines more rapidly with ra- 
dius than the composite radial surface brightness profile 
for r >0'.'27. It is therefore unlikely that the observed 
'breaks' result from the halos and structure of the ACS 
PSFs. Specifically, at z — 6 the most significant devia- 
tions in the light-profiles are seen at levels 1.5-2.0 mag 
above the Icr sky-subtraction error, and well above the 
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Fig. 7. — Composite images for Left z ~ 4, Center z~5 and 
Right 2 ~ 6 objects. Here each composite image is a stack of 30 
objects. Each stamp is 1'.'53 on a side. 




1; i 
If- I. 



^1 



Fig. 8. — Mean surface brightness profiles with a best fit Sersic 
profiles for 30 composite images at (a) z ~ 4, (b) 2 ~ 5 and (c) 
z~6, respectively. The thin dot-dash curve represents the ACS V, 
i' and z'-band PSFs, respectively, while the horizontal dashed line 
indicates the surface brightness level corresponding to the Icr sky- 
subtraction error in the HUDF images. The vertical dotted line 
marks the radius at which the profile starts to deviate significantly 
from the extrapolation of the inner r^/" profile observed at smaller 
radii. The n and ic are the best fit Sersic parameters. 

PSF wings. Each of the mean surface brightness profiles 
display a well-defined break, the radius of which appears 
to change somewhat with redshift. The vertical dotted 
lines (in Figure 6 and Figure 8) mark the radius at which 
the mean surface brightness profiles start to deviate sig- 



nificantly from the extrapolation of the 
served at smaller radii. 



r-l/" 



profile ob- 



4.3. Test of the Stacking Technique on Nearby 
Galaxies 

To test the general validity of the stacking technique 
itself on a local galaxy sample, we used surface photome- 
try from the Nearby Field Galaxy Survey (NFGS: Jansen 
etal. 2000a,b). The NFGS sample contains 196 nearby 
galaxies, that were objectively selected from the CfA red- 
shift catalog (CfA I; Davis & Peebles 1983; Huchra et al. 
1983) to span the full range in absolute B magnitude 
present in the CfA I (-14.7 <Mb< -22.7 mag). The 
absolute magnitude distribution in the NFGS sample ap- 
proximates the local galaxy luminosity function (e.g., 
Marzke etal. 1994), while the distribution over Hubble 
type follows the changing mix of morphological types as 
a function of luminosity in the local galaxy population. 
The NFGS sample (as detailed in Jansen et al. 2000a) 
minimizes biases, and yields a sample that, with very 
few caveats, is representative of the local galaxy popu- 



lation. As part of the NFGS, UBR surface photometry, 
both integrated (global) and nuclear spectrophotometry, 
as well as internal kinematics were obtained (see Jansen 
& Kannappan 2001). Here, we will concentrate on the 
fj-band surface photometry, since it is closest in wave- 
length to the rest-frame wavelengths observed at 2;~4— 6. 
Although, ideally, we would want a filter further into the 
UV, Taylor et al. (2007) and Windhorst et al. (2002) show 
that for the majority of late-type nearby galaxies, the ap- 
parent structure of galaxies does not change dramatically 
once one observes shortward of the Balmer break. Early- 
type galaxies, however, are a clear exception to this, but 
these are not believed to dominate the galaxy population 
at z~4 — 6, as discussed before. 

Figure 9 shows stacked profiles for relatively lumi- 
nous early-, spiral-, and late- type galaxies drawn from 
the NFGS. Vertical dotted lines indicate the half-light 
radii and their intersection with the profiles, the surface 
brightness at that radius. Dashed lines indicate exponen- 
tial fits to the outer portion of each profile. Figure 9 also 
shows that co-adding profiles for disparate morphologi- 
cal types and for mid-type spiral galaxies with a range 
in bulge-to-disk ratios can produce breaks in the com- 
posite profile. No such breaks are seen when the profiles 
of either early-type galaxies (E, SO) or late-type galaxies 
(Sd-Irr) are co-added. This figure shows that, if galax- 
ies at z~4 — 6 had similar morphological types as local 
galaxies, then it would be possible to produce a break in 
the profiles (as shown in Figure 6 and Figure 8), merely 
by mixing different types of galaxies. We do not believe 
that the galaxy populations at z~4 — 6 morphologically 
resemble those at low redshift. Hence, for primarily late- 
type galaxies, which dominate the faint blue galaxy pop- 
ulation at AB>24 mag (Driver etal. 1998), and which 
likely dominate the fainter end of the luminosity func- 
tion at z~4 — 6 that we sample here (Yan & Windhorst 
2004a, b), the image stacking is likely a valid exercise. 

The primary goal of this section was to show that the 
profile stacking technique is valid and can be used to 
get meaningful surface brightness profiles. We are not 
comparing our nearby sample with galaxies at z~4 — 6. 
These nearby galaxies are unlikely to be local analogues 
of high redshift galaxies. If we apply surface brightness 
dimming to UV light-profiles of these nearby galaxies, 
they would be mostly invisible to HST, and in some cases 
visible to JWST in long integration (see e.g. Windhorst 
et al. 2006) . This is another way of saying that the z ~ 4—6 
objects are truly different from z~0 objects. 

5. DISCUSSIONS 

Figure 8 shows that the mean surface brightness pro- 
files deviate significantly from an inner r^/" profile at 
radii r >0'.'27-0'.'35, depending somewhat on the redshift 
bin. These deviations appear real, with the break/point 
of departure located > 1.5-2 mag above the la sky- 
subtraction error and above the PSF-wings. In the fol- 
lowing, we discuss several possible explanations for the 
observed shapes of our composite surface brightness pro- 
files. 

5.1. Galaxies with Different Morphologies 

Our test on nearby galaxies (Figure 9) shows that, if we 
stack many galaxies with different morphologies (early- 
type, late- type or spiral galaxies), it is possible to get a 
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Fig. 9. — Stacked radial surface brightness profiles for relatively 
luminous early-, spiral- and late-type nearby galaxies drawn from 
the Nearby Field Galaxy Survey (Jansen ctal. 2000a, b). The verti- 
cal dotted line indicates the half-light radius, while the dashed line 
represents an exponential fit to the outer portions of each compos- 
ite profile. Co-adding profiles for disparate morphological types 
and for spiral galaxies with a range in bulge-to-disk ratios can pro- 
duce breaks in the composite profile. No significant breaJcs are 
seen in the outer light profiles, when the profiles of either early- 
type galaxies (E, SO) or late-type galaxies (Sd-Irr) are co-added. 

slope-change ('break') in the average surface brightness 
profile. Ravindranath et al. (2006) find that 40% of the 
brighter LBGs at 2.5 < ^; < 5 have light profiles close 
to exponential, as seen for disk galaxies, and only ~30% 
have high n, as seen in nearby spheroids. They also find a 
significant fraction (~30%) of galaxies with light profiles 
shallower than exponential, which appear to have mul- 
tiple cores or disturbed morphologies, suggestive of close 
pairs or on-going galaxy mergers. Therefore, if galaxies 
at z~4 — 6 have a variety of morphological types, then 
the shape of the average surface brightness profile that 
we see may be due to the stacking of different types of 
galaxies. Therefore, we find that the exponential and 
the flatter profiles found by Ravindranath etal. (2006) 
for galaxies at 2.5 < 2: < 5 also apply to higher redshifts 
(z>5). 

Also, we believe that it is more likely that the high red- 
shift, faint galaxy population consists primarily of small 
galaxies with late-type morphologies and with sub-L* lu- 
minosities, as seen at z~2 — 3 (Driver ctal. 1995, 1998). 
So if the z~4— 6 population consists of such a late- type 
galaxy population, then the slope-change in the light pro- 
files is likely not the result of co-adding images of objects 
with disparate morphological types. 

5.2. Central Star Formation/ Starhurst 

HST optical images of galaxies at 2; ~ 4—6 sample their 

rest-frame UV (^1200 A), where the contribution from 
the actively star-forming regions (very young, massive 
stars) dominates the UV-light. Hathi ctal. (2007) have 
shown that galaxies at 2; ~ 5-6 are high redshift starbursts 
and these galaxies have similar starburst intensity limit 
as local starbursting galaxies. Therefore, it is possible 
that galaxies at z~A—& have centrally concentrated star 
formation or starburst. This possibility is based on three 



key assumptions: (1) most of the galaxies at z~4, 5,6 
are intrinsically later- type galaxies (Driver etal. 1998; 
Steidel etal. 1999); (2) the Spectral Energy Distribution 
(SED) of these galaxies at 2 ~ 4, 5, 6 are dominated by 
early A- to late 0-type stars, respectively; and (3) there 
are no old stars with ages at 2 ~ 4 — 6 greater than 2-1 
Gyr in WMAP cosmology, respectively. 

Hunter & Elmegrecn (2006) studied azimuthally av- 
eraged surface photometry profiles for large sample of 
nearby irregular galaxies. They find some galaxies have 
double exponentials that are steeper (and bluer) in 
the inner parts compared to outer parts of the galaxy. 
Hunter & Elmegrcen (2006) discuss that this type of be- 
havior is expected in galaxies where the centrally concen- 
trated star formation or starburst steepens the surface 
brightness profiles in the center. If that is the case, then 
one might expect a better correlation between the break 
in the surface brightness profiles and changes in color pro- 
files. Unfortunately, for our sample of galaxies at z~4— 6, 
we don't have high- resolution restframe UBVcolov in- 
formation. The objects are generally too faint for Spitzer 
Space Telescope, and hence we cannot confirm or reject 
this possibility for the shape of our composite surface 
brightness profiles. 

5.3. Limits to Dynamical Ages for 2;~4, 5, 6 Objects 

The average compact 2 ~ 4 — 6 galaxy is clearly ex- 
tended with respect to the ACS PSFs (Figure 8), and 
is best fit by an exponential profile (n < 2) out to a ra- 
dius of about r ~0'.'35, 0'.'31, and 0'.'27 at 2 ~ 4, 5 and 
6, respectively. The apparent progression with redshift 
is noteworthy. The radius at which the profile starts to 
deviate from r^/" (in this case at radius r>0'.'35-0'.'27) 
may be an important constraint to the dynamical time 
scale of the system, as discussed in § 1. If this argument 
is valid, then we can estimate limits to the dynamical 
ages of 2 ~ 4, 5, 6 galaxies as follows. 

In WMAP cosmology, a radius of r>0'.'35 at z ~ 4 
corresponds to r>2.5 kpc. The dynamical time scale 
(e.g., Binney & Tremaine 1987), T^yn, goes as Tdyn = 
Cr^/'^/^/G M , where the constant C = 7r/2. For a typ- 
ical dwarf galaxy mass range of ^ 10^ — 10^ Mq inside 
r=2.5 kpc, we infer that the limits to the dynamical age 
would be Tdyn — 90-290 Myr, which is the lifespan ex- 
pected for a late-type B-star. This means that the last 
major merger that affected this surface brightness pro- 
file and that triggered its associated starburst may have 
occurred ~0.20 Gyr before z ~ 4, — assuming that the 
star-formation wasn't spontaneous, but associated with 
some accretion or a merging event. 

Table 2 shows the break-radius and inferred limits to 
dynamical ages for the z ~ 4 — 6 objects. At z ~ 5, we 
find that the limits to dynamical age at the break radius 
would be Tdyn — 70-210 Myr, which is the lifespan ex- 
pected for a mid B-star, while at z ~ 6, Tdyn — 50-150 
Myr, which is the lifespan expected for a late O early 
B-star. This means that the last major merger that af- 
fected these surface brightness profiles at z ~ 5 and 6 and 
that triggered its associated starburst may have occurred 
~0.14 and ~0.10 Gyr before z~5 and 6, respectively. 

The dynamical time is a lower limit to the actual time 
available, since it assumes matter starts from rest. Any 
angular momentum at start will increase the available 
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TABLE 2 

Dynamical Ages for z^i — Q objects in the HUDF 



Redshift 
z 


"Break" Radius" 
(arcsec) 


"Break" Radius'' 
(kpc) 


Dynamical Age'^ 
iXdyn ) 


4 


0.35 


2.5 


0.09-0.29 Gyr 


5 


0.31 


2.0 


0.07-0.21 Gyr 


6 


0.27 


1.6 


0.05-0.15 Gyr 



^ From composite surface brightness profiles (Figure 6 and Figure 8). 
^ Radius in kpc corresponding to radius in arcsec at given redshift. 
^ If "break radius" interpreted as indicator of dynamical age. 



time. The best-fit SED age from the GOODS HST and 
Spitzer photometry on some of the brighter of these ob- 
jects — using Bruzual & Chariot (2003) templates — is 
in the range of about ^150 650 Myr (Yan ctal. 2005; 
Eyles etal. 2005, 2007), the lower end of which is con- 
sistent with our limits to their dynamical age estimates, 
while the somewhat larger SED ages could also be af- 
fected by the onset of the AGB in the stellar population 
increasing the observed Spitzer fluxes and hence possibly 
overestimating ages (Maraston 2005). Our age estimates 
for 2^4 — 6 are consistent with the trend of SED ages 
suggested for 2~7 (Labbe etal. 2006). It is noteworthy 
that, given the uncertainties, the two independent age 
estimates are consistent. If our limits to dynamical age 
estimates for the image stacks are thus valid, they are 
consistent with the SED ages, and point to a consistent 
young age for these objects. 

Furthermore, the presence of young, massive late O- 
early B-stars at z ~ 6 has implications for the rcionization 
of the universe. From observations of the appearance of 
complete Gunn-Pctcrson troughs in the spectra of z> 5.8 
quasars (Fan etal. 2006), we know that the epoch of 
reionization had ended by 2; ~ 6. From the steep (a=-l .8) 
faint-end slope of the luminosity function of z ~ 6 galax- 
ies, Yan & Windhorst (2004a;b) concluded that dwarf 
galaxies, and not quasars, likely finished rcionization by 
2 ~ 6. Should the present interpretation of their light 
profiles be correct, then it would appear to add support 
to this picture, in the sense that such objects are dom- 
inated by B-stars and did not start their most recent 
major starburst long before 2;~6. 

6. SUMMARY 

We used the stacked HUDF images to analyze the av- 
erage surface brightness profiles of z ~ 4 — 6 galaxies. 
Our analysis shows that even the faintest galaxies at 
2; ~ 4 — 6 are resolved. This may have implications on 
the stellar density and its relation to the stellar density 
in present-day galaxies. We also find that the average 
surface brightness profiles display breaks at a radius that 
progresses toward lower redshift from r ~0'.'27 (1.6 kpc) 
at 2~6 to r :i0'.'35 (2.5 kpc) at z~4. 



The shape of the radial surface brightness profile that 
we observe could result from a mixture of different mor- 
phological types of galaxies, if they exist at z ~ 4 — 6, 
because we can produce similar breaks in the surface 
brightness profiles when we mix different types of nearby 
galaxies. Alternatively, if these galaxies are dominated 
by a central starburst then they could show such dou- 
ble exponential-type profiles, as discussed by Hunter & 
Elmegreen (2006). In a third scenario, if the galaxies 
at 2; ~ 4 — 6 are truly young and mostly late-type, the 
outer profiles seen in our mean radial surface brightness 
profiles at z ~ 4 — 6 bear the imprint of the hierarchical 
build-up process, and are still dominated by infalling ma- 
terial, which is not detectable in the individual HUDF 
images of these faint objects. We have estimated limits 
to dynamical ages from the break radius at z ~ 4, 5, 6, 
very roughly as ~0.20, 0.14 and 0.10 Gyr, respectively, 
and those ages are similar to the SED ages inferred at 
z~4-6 (Yan et al. 2005; Eyles et al. 2005, 2007), and con- 
sistent with SED ages suggested for ^ ~ 7 (Labbe et al. 
2006). Hence, at z~4, 5,6, the last major merger that 
affected the surface brightness profiles that wc observe, 
and that triggered the observed star- burst, may have oc- 
curred respectively ~0.20, 0.14 and 0.10 Gyr earlier, or 
very approximately at z ~ 4.5,5.5,6.5. This would be 
consistent with the hierarchical assembly of galaxies and 
with the end of reionization, since it would imply that 
from z~4 to z~6, the SEDs become progressively more 
dominated by late-B-late-0 stars. This implies that the 
sub-L* (i.e. dwarf) galaxies may have produced sufficient 
numbers of energetic UV photons to complete the reion- 
ization process by z~6, as Yan & Windhorst (2004a, b) 
suggested. It will be imperative to study with future in- 
struments fike HST and JWST (Windhorst etal. 
2006, 2007) whether the dominant stellar population in- 
deed changes from late-O-early-B at z ~ 6 (i.e. capable 
of reionizing) to mid- to latc-B at z ~ 4 — 5 (i.e. capa- 
ble of maintaining reionization), and to what extent the 
intrinsic sizes of these faint objects will ultimately limit 
deep JWST surveys. 
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